function y = rst_h_sta(temp, rym)
    Fs = 8192;
    wavef = [1046, 1174, 1318, 1396, 1566, 1760, 1976];
    x = linspace(0, 2 * pi * rym, floor(Fs * rym));
    y = sin(wavef(temp) * x);
end

